Imaging correlated wave functions of few-electron quantum dots: 
Theory and scanning tunneling spectroscopy experiments 
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We show both theoretically and experimentally that scanning tunneling spectroscopy (STS) im- 
ages of semiconductor quantum dots may display clear signatures of electron-electron correlation. 
We apply many-body tunneling theory to a realistic model which fully takes into account correlation 
effects and dot anisotropy. Comparing measured STS images of freestanding InAs quantum dots 
with those calculated by the full configuration interaction method, we explain the wave function 
sequence in terms of images of one- and two-electron states. The STS map corresponding to double 
charging is significantly distorted by electron correlation with respect to the non-interacting case. 
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I. INTRODUCTION 

Scanning tunneling microscopy (STM) and spec- 
troscopy (STS) is a key tool in nanoscience, allowing for 
both manipulation of nano-objects and access to their 
energy spectrum and wave function (WF). The class of 
systems under study is broad and loosely defined, ranging 
from nanotubes to quantum dots (QDs) and molecules. 
Their common features are the possibility of achieving 
relatively good electrical insulation, the discreteness of 
their energy spectrum, and the manifestation of Coulomb 
blockade phenomena at low temperatures. In these few- 
body systems electron-electron interaction may play a 
major role, as it is immediately apparent e.g. from single- 
electron charging experiments. 1 

Some of us recently suggested that the WF imaging 
technique of STS cou ld be a sensitive and direct probe 
of electron correlation! 2 * 3 * 4 * 5 ^ Despite recent experimental 
evidence that few-body semiconductor QDs are strongly 
affected by electron correlation, as seen both in inelastic 
light scattering^ and high source-drain^ spect rosc opies, 
all QD WF images obtained so far, both in rea l 8 * 9 * 10 * and 
reciproca l 11 * 12 **^ space, were basically interpreted in 



terms of independent-electron WFs. In this paper we 
apply many-body tunneling theory to a realistic model 
which fully takes into account the combined effect of 
dot anisotropy and Coulomb interaction and predict that 
STS maps can be strongly distorted by correlation ef- 
fects. These calculations are validated by a recent low- 
temperature-STS experiment performed by a few of us 
on few-electron InAs QDs at the University of Hamburg, 
which was reported in Ref. [4] and that we here review 
providing further details and discussion of the theoreti- 
cal model. 

While in this paper we focus exclusively on semicon- 



ductor QDs, we stress that the main ideas regarding how 
electron correlation affects STS could in principle be ap- 
plied also to short carbon nanotubes as well as single 
molecules. 

The plan of the paper is as follows: After a re- 
view of conceptual frameworks used to understand STS 
(Sec. |ll|, we illustrate our full configuration int erac tion 
(FCI) method of calculation of STS images (Sec. Ill) and 
theoretically focus on the double charging of a QD in var- 
ious realistic regimes (Sec. |FV| . We discuss the evidence 
of correlation effects provided by the Hamburg experi- 
ment (Sec.[V]), and eventually reconsider our findings in 
the Conclusions (Sec. VI). 



II. 



THEORY OF THE STS IN THE PRESENCE 
OF ELECTRON CORRELATION 



In order to understand the origin of correlation ef- 
fects in STS, it is useful to critically reconsider the basic 
picture of STS performed in the presence of Coulomb 
blockade. Figure 1 displays schematically the typical 
energy landscape along the tunneling direction seen by 
an electron flowing from the STM tip to a drain lead 
(backgate in Fig. 1) through the QD. The QD is well 
isolated from tip and backgate by two tunneling barri- 
ers (the vacuum and an insulating layer, respectively), 
so the width of its discrete energy levels is negligible. 
At low temperature transport occurs mainly due to res- 
onant tunneling through the QD chemical potentials 
/jl(N) = Eq(N) — Eq(N — 1), which are depicted as seg- 
ments in the QD region of Fig. 1 [N is the number of 
electrons filling in the QD at equilibrium and Eq(N) is 
its ground-state energy] P^* If the tip-backgate voltage V 
is small, current may flow or not (Coulomb blockade) de- 
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FIG. 1: Simplified energy landscape along the direction of 
tunneling through a quantum dot for a typical scanning tun- 
neling spectroscopy measurement. Ey is the Fermi energy 
of the STM tip. Upper right inset: Corresponding idealized 
dl/dV vs. V plot (ignoring all broadening mechanisms but 
resonant tunneling). 



used in the iV- dep endent computation of self-consistent 
SP orbitals Vv 19 (ii) The interacting TV-electron ground 
state \^n) is m general a linear superposition of many 
electronic configurations (Slater determinants), l^f^), 
with expansion coefficients cf , \$n) = J2i c i* 
while mean-field theory approximates \^n) with a single 
Slater determinant, {cf = Sn): this approximation 

turns out to be poor in strongly correlated regimes.^ 

In order to circumvent the above difficulties we recall 
from many-body tunneling theory that the differential 
conductance is p roportional to the interacting local den- 
sity of states p2 



dI/dV(E F ,r) oc -— Im0(r, r; E F ), 



(i) 



where Q(r,r;E F ) is the energy- and space-resolved ex- 
act zero-temperature retarded Green's function (or one- 
electron propagator).^ The quantity on the r.h.s. of 
Eq. ([!]) (also known as spectral density) may be regarded 
as the squared modulus of a quasi-particle WF: 2 



pending if some value of N exists such that fi(N) enters 
the transport energy window. By increasing V one is able 
to widen the transport window, causing a step in the cur- 
rent I to occurr (a peak in the differential conductance 
dl/dV in the inset of Fig. 1) each time the chemical po- 
tential /i(iV) for a new electron number N' enters the 
window. 

In the standard STS mean- field theory, 16 the chemi- 
cal potentials /i(iV)'s are simply Hartree-Fock or Kohn- 
Sham single-particle (SP) self-consistent levels, split in 
energy by Coulomb and exchange interactions in addi- 
tion to the effect of quantum confinement. The energy- 
and space-resolved local density of SP states, n(e,r), 
is just given by a sum over the occupied SP orbitals 
^a(r), n(e,r) = £) a |^ a (r)| 2 S(e - e a ), where e a is the 
energy of the ath SP level. A general result of many- 
body time-dependent perturbation theory, valid to first 
order in the tunneling matrix element,^ is that dl/dV 
is proportional to the local electron density. Therefore, 
in the framework of mean-field theory, the local value of 
the differential conductance, dl j dV (E F , r) , is given by a 
sum over the square moduli of SP orbitals whose energies 
fall in the transport window Se around the Fermi energy 
E F , dI/dV(E F ,r) oc ^2 Se \ifj a (r)\ 2 . If the energy resolu- 
tion is sufficient, the space-resolved map of dI/dV(E F , r) 
provides the image of the square modulus of a single 
SP orbital, namely the highest-energy occupied QD or- 
bital, akin to the HOMO in molecules. Besides, the STS 
energy-scan mode exhausts the Aufbau filling sequence of 
the lowest-energy SP levels of the dot. 

Since in general the QD energy spectrum strongly de- 
pends on electron-electron interaction and consequently 
on TV, 18 the above mean- field theory suffers a few serious 
drawbacks: (i) Each tunneling event involves the transi- 
tion between QD ground states with TV and N + 1 elec- 
trons, therefore it is unclear which value of N should be 



lmg(r,r;E F ). 



(2) 



The quasi-particle WF (^qd(^) is the natural generaliza- 
tion of the SP WF to strongly correlated regimes. In 
the non-interacting limit ^qd(^) —> as we show 

below. In practice, the calculation of (^qd(^) requires 
the knowledge of the (configuration interaction) expan- 
sion coefficients and cf~ x of the interacting ground- 
states |\£jv) and I^jv-i) with both N and TV — 1 elec- 
trons, respectively, according to the formula ^qd(^) = 
. cf r_1 *cj vr '0 Q; (^j)(r), where a(i,j) is a SP quantum 
index depending on both Slater determinants l^^ -1 ) 
and |$^)- Distortion effects of (/?qd(t) with respect to 
ip(r) are due to the quantum interference among differ- 
ent SP orbitals ip a (ij)( r )' I n the non-interacting limit, 
both N- and (N — l)-electron ground states are single 
Slater determinants, cf~ x = Sn, cj^ = Sji, and therefore 
^qd(^) reduces to the simple SP orbital ^ a (i,i)( r )- 



III. PREDICTING STS MAPS FROM FULL 
CONFIGURATION INTERACTION 

We consider a two-dimensional QD with parabolic 
lateral confinement, which is a commonly accepted 
approximation^ in the framework of the envelope- 
function description of SP states: 



H (i) = 



Pi 



2 2\ 



(3) 



2m* 2 

In the SP Hamiltonian for the zth electron of Eq. ([3| 
the lateral confinement is different in the x and y direc- 
tions (with corresponding confinement frquencies uo x and 
ujy, respectively): Such an elliptical confinement mimics 
the combined effects of geometrical deviations from per- 
fect circularity and/or atomistic effects 23 due to strain, 
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piezolectric fields, interface matching, which lower the 
symmetry point-group from D^h (circular case) to C2 V 
[the symmetry point-group of the Hamiltonian ([3]), D^h^ 
is actually slightly larger than the C2 V group]. In addi- 
tion, we include into our model the full Coulomb inter- 
action among the N electrons populating the dot: 

N i 2 

In Eqs. (|3][4]) e and k are respectively the electron charge 
and static relative dielectric constant of the host semi- 
conductor, p is the canonically conjugated momentum of 
position r = (x,y), m* is the electron effective mass. 

We solve numerically the few-body problem of Eq. Q 
for the ground states of N and N — 1 electrons by means 
of the FCI method, that we successfully applied in pre- 
dicting QD transport^ and RamaiP spectra (for full de- 
tails on our FCI method, its performances, and ranges 
of applicability, see Ref . l20|) . Briefly, we expand the N- 
electron interacting ground state \^n) on the basis of 
the Slater determinants obtained by filling in with 

N electrons in all possible ways a subset of the SP or- 
bitals, eigenstates of Hamiltonian (13b - On the Slater- 
determinant basis, the Hamiltonian Ml) is a large, sparse 
matrix that we diagonalize by means of a parallel state- 
of-the-art Lanczos code. The expansion coefficients cf , 
which are the output of the FCI calculation, are eventu- 
ally used to build the quasi-particle WF (^qd(^)- 

IV. QUASI-PARTICLE IMAGES IN REALISTIC 
QUANTUM DOTS: THEORY 

We here consider quasi-particle WF images for realistic 
QDs predicted by the theory of Sees. |H]|III[ focusing on 
the simplest case where electron correlation becomes rel- 
evant, namely the tunneling transition N = 1 — > N — 2. 
In fact, the image of the first charging process N = — > 
N = 1 is simply given by the lowest-energy SP orbital. 

We first investigate the case of a circular QD by moni- 
toring the evolution of the square modulus of ^qd(^) as 
a function of the lateral confinement energy Hluo = fow x — 
hujy. On a general basis, we expect that correlation ef- 
fects are negligible in the non-interacting limit Hluq — > oo 
while they dominate in the opposite limit Hljo — > 0. In 
fact, the SP [Eq. and the Coulomb [Eq. Q] terms 
of the interacting Hamiltonian scale differently with hwo , 
the former as ~ Jtluq and the latter as ~ V^^o (Ref. [2j). 
When ftwo — > oo the SP term dominates with respect to 
the Coulomb term, and the ground state is essentially 
a single Slater determinant, while interaction effects can 
be regarded as a perturbation: in this limit mean-field 
theory is expected to give correct predicitons. In the op- 
posite limit, instead, the Coulomb term very effectively 
mixes many Slater determinants, and significant distor- 
tions of ^qd(^) with respect to the non-interacting case 
are expected. 




FIG. 2: Calculated STS maps for the ground state — > ground 
state transition N = 1 ^ N = 2, as a function of lateral con- 
finement energy hu>o and dot anisotropy. (a-c) Dependence of 
the STS image on the lateral confinement energy Hujq in the 
case of a circular dot. Plots (a-c) correspond to Huoq = +oo, 
0.5, 0.01 meV, respectively, with ujq — uj x = oj v . (d-e) De- 
pendence of the STS image on the dot anisotropy. Plots (d-e) 
correspond to (jj x /u y = 1, 1.1, 1.5, respectively, with fwj y = 
2 meV. The lateral extension of all plots is 4 x 4 units of 
the characteristic lateral extension of the harmonic oscillator 
^qd = [h/ (m* Uy)] 1 / 2 , and GaAs parameter are used through- 
out the paper. Plot heights are renormalized arbitrarily. The 
absolute norms of STS maps are given in Table |I1 



Figures 2 (a-c) display the maps of |(/?qd(^)| 2 in the 
(x,y) plane for decreasing values of Hluq [from top to 
bottom, hujQ — > +oo (a), Houq = 0.5 meV (b), Huuq = 0.01 
meV (c)]. In the non-interacting case [Fig. 2(a)], the STS 
map is just a replica of the SP Is orbital, eigenstate of the 
2D harmonic oscillator (a 2D gaussian), corresponding to 
the injection of a second electron into the same Is level 
occupied by the first electron with opposite spin, accord- 
ing to Pauli exclusion principle. As Hujq is decreased, we 
see that a significant portion of weight is moved from the 
image center into an outer ring [Fig. 2(b)]. In a very shal- 
low dot [Fig. 2(c)], the dot center is completely emptied 
and the STS image now looks like a donut. The latter 
case of strong interaction shows that the STS image can 
be very different from that expected by a naive applica- 
tion of the Aufbau principle, due to the hybridization of 
higher-energy SP orbitals with the Is WF. 

In Fig. 2 the length unit is the characteristic lateral ex- 
tension of the harmonic oscillator, ^qd = [^/(ra*^)] 1 / 2 , 
which depends on ftw y , and ^qd(^) is arbitrarily normal- 
ized. Another effect of correlation is the loss of absolute 
weight of the quasi-particle WF as hwo is decreased. 2 
Table [i] shows the norm of ^qd(^) for the cases dis- 
played in Fig. 2. We see that the reduction of the con- 
finement energy Hujq in a circular QD is associated to 
a dramatic weight loss. This trend may be seen as a 
signature of the Wigner crystallization of the (N — 1)- 
and A/"-electron ground states which manifests itself as 
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TABLE I: Norm of the quasi-particle, f dr |^qd(V)| 2 , for 
different ground-state — > ground-state tunneling transitions 
N — 1 — > A/" = 2. The weight ranges between and 1. 



Case a 


Quasi-particle weight 


Case 


Quasi-particle weight 


(a) 


1.00 


(d) 


0.794 


(b) 


0.633 


(e) 


0.793 


(c) 


0.109 


(f) 


0.720 



a The case refers to the calculated STS image displayed in Fig. [2] 



an increased "rigidity" of the states opposing electron 
injection. 

We now switch to consider the effect of dot anisotropy, 
since the dependence of STS images on dot ellipticity 
was found to be a major issue in the Hamburg experi- 
ment (cf. Sec. |V| below) Pin Fig. 2(d) we start showing 
the predicted STS map for a circular QD with Hujq = 
2 meV, and then modify its ellipticity by increasing the 
lateral frequency ratio cu x /cu y , going from co x /co y = 1.1 
[Fig. 2(e)] up to u) x /u y = 1.5 [Fig. 2(f)]. By doing so, 
as we move downwards in the right column of Fig. 2, we 
expect to "squeeze" the quasi-particle WF along the x di- 
rection, which is indeed observed for both single and dou- 
ble charging processes. However, in the present double 
charging case, we also see an unexpected effect, i.e. the 
quasi-particle WF develops two peaks along the y-axis. 
This surprising distortion is due to the destructive inter- 
ference between Is and Id states of the harmonic oscil- 
lator along y (belonging to the same representation A of 
the D2h group), which is a correlation effect. This can 
also be seen as a manifestation of the general statement 
that the importance of correlation increases as the system 
dimensionality is reduced (from 2D to ID). 

To conclude this section, we mention that another is- 
sue relevant to STS images is the type of dielectric en- 
vironment felt by QD electrons.^ Specifically, the dielec- 
tric mismatch between vacuum and InAs in the Ham- 
burg experiment affects both the SP confinement poten- 
tial (self-polarization effect) and the electron-electron in- 
teraction (interaction with surface image charges of like 
sign), therefore changing the relative importance of their 
effect on the ground state. Qualitatively, we expect that 
changing the dielectric environment causes effects similar 
to those obtained by modifying Huoq. 



V. STS IMAGING OF FEW-ELECTRON 
MBE- GROWN QUANTUM DOTS: EXPERIMENT 

We studied strain-induced InAs QDs grown on n- 
doped GaAs(OOl) substrates by molecular beam epitaxy 
(MBE) An undoped tunneling barrier 5 nm thick insu- 
lates the freestanding QDs from a n-doped GaAs buffer 
layer 200 nm thick, acting as a backgate for the current 
flowing from the STM tip through the QDs. The samples 
were transferred from the MBE into the STM chamber 




FIG. 3: dI/dV(E,r) spatial maps recorded at the voltages 
840 (map a), 1040 (b), 1140 (c), and 1350 (d) mV, respec- 
tively, for a representative QD. The size of all maps is 30 x 
30 nm. White (black) colour stands for high (low) values of 
dl/dV. /stab = 100 pA, V stah = 1500 mV, V mod = 4 mV. 



in ultra-high vacuum without being exposed to the air by 
means of a mobile system, in ordered to avoid contami- 
nation, and STM was operated at a background pressure 
smaller than 10 -10 mbar and at T = 6 K with maximum 
energy resolution of Se = 2 meV. 25 STM images were 
taken in constant-current mode, with a typical sample 
bias in the range of 2-4 V and a tunneling current of 20- 
40 pA. A lock-in technique (modulation voltage V mo a in 
the range of 5-20 mV) was used to record dI/dV(E,r) 
and WF mapping was carried out over a specified area 
by stabilizing the tip-surface distance in each point r at 
voltage F s tab and current / s tab, switching off the feed- 
back and recording a dl/dV curve from V^tart to V en( ± 
(K$t art < Ktab)- As a result, WF mapping produces a 
3D array of dl/dV data, which allows obtaining spatially 
resolved dI/dV(E, r) images at different values of sample 
voltages. 

In Fig. 3 we show dI/dV(E,r) maps of a representa- 
tive QD measured at fixed voltages corresponding to four 
clearly localized QD WFs, labeled (a), (b), (c), and (d), 
respectively, in order of increasing energy!^ We observed 
a pronounced shape anisotropy for all dots, which have 
a pyramidal shape wi th we ll-defined facets, consistently 
with previous findings! 10 * 26 * This is evident also from the 
elongation of images (b), (c), and (d) of Fig. 3 along 
the [110] direction. In detail, WF maps of Fig. 3 show 
the following approximate symmetries, going from low 
to high-energy: s-like for (a), p-\ike for (b) and (c), and 
possibly (d). The (d) image is somehow blurred, likely 
because states of the wetting layer overlapping in energy 
with (d) significantly contribute to the spectral density. 
See Ref. |4] for further discussion. 

As expected for the square moduli of the two lowest- 
energy SP Is and 2p orbitals, images (a) and (b) exhibit 
a roughly circular symmetric intensity distribution and 
elongation along the [110] direction with a node in the 
center, respectively (Fig. 3). Unexpectedly, image (c) 
shows again a p-\ike symmetry in the [110] direction, as 
before, instead of [110] as expected for the second 2p 
orbital. 10 As a consequence, it is not possible to explain 
the WF sequence [and map (c) in particular] in terms of 
SP orbitals, since in this case we would expect the ap- 
pearance of either a single 2p or two 2p states elongated 
in the [110] and [110] directions, respectively. On the 
other hand, from a mean-field point of view the charg- 
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ing of the same 2p orbital with a second electron can be 
excluded since a replica is not observed for the Is orbital. 

On this basis, we believe that the understanding of 
the data of Fig. 3 must rely on the theory developed in 



Sec. [TV] taking into account the combined effect of elec- 
tron correlation, dot anisotropy, and dielectric mismatch. 
Specifically, we assign maps (a) and (b) to the tunneling 
events N = —> N = 1 corresponding to resonances 
through Is and 2p SP states, while we associate image 
(c) of Fig. 3 to the ground state — > ground state charg- 
ing process TV = 1 — > N = 2, showing features similar to 
those predicted in Figs. 2(e-f). Further support for this 
interpretation is reported in Ref. |4j 



VI. CONCLUSIONS 

In this paper we have focused on QDs where electron- 
electron interaction may be relevant, showing that the 
STS technique is sensitive to quasi-particle WFs and that 
their images can be greatly affected by electron corre- 



lation. On the basis of FCI calculations of STS maps 
which fully take into account correlation effects and dot 
anisotropy, we have been able to understand measured 
WF images of freestanding self- assembled InAs quantum 
dots. We identified ground- and excited-state WFs corre- 
sponding to the injection of a first and a second electron 
into the QD. Correlation effects are found to distort the 
WF corresponding to double charging in an essential way. 
We believe our findings may be relevant to a broader class 
of experiments, including STS on carbon nanotubes and 
single molecules. 
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